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We study theoretically the electrical current and low-frequency noise for a linear Josephson junc- 
tion structure on a topological insulator, in which the superconductor forms a closed ring and 
currents are injected from normal regions inside and outside the ring. We find that this geometry 
offers a signature for the presence of gapless ID Majorana fermion modes that are predicted in the 
channel when the phase difference ip, controlled by the magnetic flux through the ring, is ir. We 
show that for low temperature the linear conductance jumps when ip passes through ir, accompanied 
by non-local correlations between the currents from the inside and outside of the ring. We compute 
the dependence of these features on temperature, voltage and linear dimensions, and discuss the 
implications for experiments. 
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I. INTRODUCTION 

There is presently a major effort in condensed mat- 
ter physics to demonstrate the unique properties of 
Majorana fermion quasiparticle states associated with 
topological superconductivity.^ A promising approach 
is to utilize proximity effect devices that combine or- 
dinary superconductors with topological insulators or 
other strong spin-orbit materials to achieve topological 
superconductivityPHU Recent experiments on semicon- 
ductor nanowires coupled to superconductors observed 
zero-energy tunneling resonances that have been inter- 
preted as Majorana bound statesPHH] The original pro- 
posal involved a superconductor coupled to the surface 
of a topological insulator (TI). It was shown that a vor- 
tex in the superconductor is associated with a Majorana 
bound state, and that a linear Josephson junction (JJ) 
exhibits gapless ID Majorana fermions when the phase 
difference is ip = w™ Supercurrents in TI JJ devices have 
recently been observed.^HUl Unusual behavior, including 
a smaller than expected critical current normal resistance 
product, as well as an anomalous Fraunhofer diffraction 
pattern, has been interpreted as evidence for ID Majo- 
rana fermions along the channel between the supercon- 
ductors.^ However, the connection between these obser- 
vations and Majorana fermions is indirect. 

A difficulty with critical current measurements is that 
the predicted Majorana behavior is only manifest when 
the phase difference ip is close to n. In order to isolate the 
properties of the gapless Majorana mode, it is necessary 
to control the phase. Moreover, the supercurrent car- 
ried by the junction includes contributions from gapped 
states. The gapless mode leads to only a weak singularity 
in the current phase relation at ip — ir. Another possible 
experiment would be to use a ring geometry and to tun- 
nel into the junction region from another contact. Similar 
experiments on ID SNS junctions revealed the expected 
collapse of the minigap in the normal region for ip = 7rP^ 
Such experiments on a TI JJ device could demonstrate 
the closing of the gap, but they would not distinguish an 
even and odd number of gapless channels. The unique 




FIG. 1: (a) Three-terminal ring geometry, in which the phase 
difference across a linear S-TI-S Josephson junction is con- 
trolled by the magnetic flux <3> through the ring. Measurement 
of the current and low-frequency noise provides a signature 
sensitive to Majorana fermion modes on the junction. 

feature of the ID Majorana mode is that as the phase is 
tuned through ir, the system must pass through a state 
where the Majorana mode is transmitted perfectly along 
the channel - even in the presence of strong disorder. This 
leads to a quantized thermal conductance that in princi- 
ple probes the central charge c = 1/2 associated with the 
Majorana mode, but would be difficult to measure. 

In this paper we show that measurement of current 
and low-frequency current noise injected from the nor- 
mal surface states at the ends of a TI JJ give rise to a 
unique signature that is sensitive to the perfectly trans- 
mitted gapless Majorana mode at ip = ir. Our setup, 
shown in Fig. [I] consists of a superconducting ring on 
the surface of a topological insulator with a linear chan- 
nel. The phase difference across the junction is controlled 
by the magnetic flux through the ring, ip = 2ir^/(f>Q, with 
cf>o = ft./2e. In addition, we introduce electrical contacts 
to the normal TI surfaces inside and outside the ring, as 
well as to the superconductor. This three-terminal ge- 
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ometry allows the average currents (I a ) (a = 1,2) and 
the correlations (5I a SI b ) to be measured as a function of 
the voltages V a relative to the grounded superconductor. 
We find that for a long junction at low temperature the 
linear conductance exhibits sharp steps when ip passes 
through odd multiples of tt. This singular behavior is 
a direct consequence of the gapless Majorana mode at 
ip = 7T. When the electrical contacts consist of only a 
single transmitting channel, the magnitude of these steps 
is 2e 2 /h. In the more realistic situation where there are 
multiple channels in the leads the singular behavior is 
still present, though the size of the step is reduced. 

In addition, we find signatures for the Majorana mode 
in the diagonal and cross correlations of the current noise. 
The cross correlation (SI^-SI 2 ) exhibits a narrow peak 
when ip ~ tt due to the Majorana mode. In the single 
channel limit, the height of the peak in the zero tem- 
perature limit is universal. For many channels, the peak 
height is suppressed. The magnitude and sign of the 
peak, however, is predicted to be related to the size of the 
steps measured in the average current. Finally, we pre- 
dict the diagonal noise correlation (6 1 1 6 1 1 ) also exhibits 
a sharp peak at tp = n. Unlike the singular behavior of 
the average conductance and the cross correlation, how- 
ever, the size of the peak is not suppressed when there 
are many channels, and provides a more robust signature 
of the Majorana mode. 

Our results are related to earlier work concern- 
ing resonant transmission through OD Majorana bound 
statesP^K2llln particular, Law, Lee and Ng^B showed that 
tunneling into a Majorana bound state at T = leads to 
a zero-bias resonant conductance 2e 2 /h associated with 
perfect Andreev reflection of a single channel. In our ge- 
ometry there are no discrete Majorana zero modes. How- 
ever, when the magnetic flux, rounded to the nearest in- 
teger multiple of 4>o is odd, there is effectively a Majorana 
zero mode inside the ring, but it is strongly coupled to 
the continuum of states in the surface region coupled to 
the lead. In this regard, the linear JJ exhibits a behavior 
similar to the topological transition in a ID topological 
superconductor^. In our case, the transition between 
the topological and non-topological phase, which occurs 
at tp — tt is controlled by the magnetic flux. Noise cor- 
relations associated with tunneling into Majorana bound 
states have also been studied, and the noise correlations 
for ip — tt resembles the noise correlations that have been 
studied for coupling to chiral Majorana fermion modes 
associated wit h mag net-superconductor interfaces on the 
surface of a TI P 1 * 28 ! An important difference between the 
present work and these earlier works, however, is that in 
our geometry, the superconducting phase ip, controlled 
by the magnetic flux through the ring, provides an ac- 
cessible knob for controlling the coupling between the 
counterpropagating chiral Majorana modes. This makes 
it possible to tune through the quantum critical point 
at ip = 7r that separates topological and non-topological 
phases. An advantage of the present approach is that 
the singular behavior of the current and noise near tp = tt 



provides a distinctive signature for the Majorana physics. 

The organization of our paper is as follows. We first 
present in section [TT| the specifics of our model for the 
Josephson junction, focusing on a simple tunneling prob- 
lem for the Majorana channel and properly treating the 
interchannel reflection of the remaining modes in the 
leads. Section |III| focuses on experimental signatures 
of the gapless Majorana channel: first we discuss the 
ideal case of single-channel leads and then we expand 
out to the more realistic multi-channel case, extracting 
the terms in the multi-channel current and noise which 
show singular behavior as tp winds through tt. Finally, 
we conclude with a discussion of experimental parame- 
ters and feasibility. Additional appendixes [A| through [C| 
provide complete derivations of the multi-channel observ- 
ables and their single- and many-channel limits. 

While this manuscript was in the final stages of prepa- 
ration we received a preprint by Diez, et al. that presents 
an analysis of a related Josephson junction geometry for 
topological superconductorsP3 



II. MODEL SYSTEM 

In Ref. [3] it was shown that a Josephson junction on 
the surface of a topological insulator exhibits a gapless 
Majorana mode when the phase difference is equal to tt. 
At that point, there is a single, one-dimensional Majo- 
rana mode that is transmitted perfectly along the length 
of the junction. For ip different from tt, a gap opens in 
that mode. In the following we consider the low energy 
limit for phase difference close to tt, so that the trans- 
mission across the junction is dominated by the single 
gapless Majorana mode, which can be described by a 
simple one-dimensional model. As indicated in Fig. [2] 
this mode couples to a single Majorana mode in the con- 
tacts inside and outside the ring. In general, the contacts 
will involve many additional channels, so it is necessary 
to introduce a general scattering matrix that relates the 
incident channels and the transmitted channel. 

We will begin with a discussion of the one-dimensional 
transmission problem in |II A| and in |II B| we introduce 
the general scattering matrix, from which we compute 
the current and noise as functions of the phase difference 
ip, temperature T and the voltages on the inner and 
outer contacts relative to the grounded superconductor. 



A. One-Dimensional Model 

In Ref. [3] a TI JJ was described by modeling the 
TI surface state by a single massless 2D Dirac fermion 
coupled to a superconducting pairing potential A(y) = 
Ao eicp(isgn(y)ip/2). For tp ~ tt there are quasiparti- 
cle states with E ~ bound to the interface y = 
that are described by a two band Hamiltonian, H = 
(7 L 7 l R )7~L( r y L , r y R ) T , where the one-body Bogoliubov de 
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FIG. 2: Schematic of the ID Majorana modes that propa- 
gate along the junction for tp ~ tt, connecting reservoirs on 
the inside and outside of the ring. A mass term Ao cos(<p/2) 
couples the counterpropagating Majorana modes. The trans- 
mitted mode at tp = -n defines a single Majorana mode in both 
contacts. The other Majorana modes in the contacts will be 
reflected, and characterized by reflection matrices R 1 and R? . 



Gennes Hamiltonian 



Ti. = —ihvpa z d x + m(x)a v . 



(2.1) 



Here 7° (x) = 7°(x)t are ID Majorana fermion operators 
with a = L, R. In this Majorana basis the one body 
Hamiltonian % exhibits particle-hole symmetry {H, EE} — 
with 3 = K, complex conjugation. The mass term is 
m = AoCOS^/2. Importantly, m changes sign when tp 
advances by 2-7T. In our ring geometry, this means that 
sgn(m) = (— 1)", where v — <&/4>o, rounded to the nearest 
integer. The mass term violates time reversal symmetry, 
expressed by [H,Q] = with = ia v K. For tp. = tt, 
m = 0, there are uncoupled counterpropagating chiral 
Majorana fermion modes on the ID interface. To model 
the ends of the junction, we suppose that A(x,y) varies 
adiabatically as a function of x and smoothly goes to zero 
in the lead regions with \x\ > L/2. In this case the ID 
Majorana modes on the channel evolve into one of the 
many propagating channels in the leads. In the spirit 
of the Landauer-Biittiker approach,^ we focus on this 
single channel and arrive at a ID model for the finite 
length J J coupled to the leads described by Eq. (2.1 1. 



This defines a scattering problem for the chiral Majorana 
fermion modes incident from the leads. 

This scattering problem can be characterized by a 2 x 2 
S matrix, 



(2.2) 



where t E ^ and r 1 ^ describe the amplitudes for trans- 
mission and reflection of quasiparticles with energy E 
incident from the left (right) side. Se obeys a number of 
general constraints. Unitarity requires \r E \ 2 + |ij|| 2 = 1 
and \t E \ 2 — \t 2 E \ 2 = \t E \ 2 . Particle-hole symmetry re- 
quires S-e = S E . The scattering problem is easily solved 
for the simple model m{x) = 8(L/2 — |x|)Ao cos <fi/2. 
This model has a mirror symmetry (x — > —x), under 



-r\ = Te- We find 
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Tivfh cosh kL — iE sinh kL ' 

m sinh kL 
TivpK cosh kL — iE sinh kL ' 



(2.3) 
(2.4) 



where k = \JE 2 — m 2 / (Tivp). At E — 0, Sq is real, and 
is characterized by 
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1/ cosh(m/Ae) 
tanh(m/Ae), 



(2.5) 
(2.6) 



where we have defined Ae = Tivf/L. For m ^S> Ae, 
ro = ±1. Importantly, when m changes sign, tq changes 
sign and must pass through zero, at which point the 
transmission at E = is perfect, i.e., \t \ = 1. This 
property is more general than our specific model. It is 
related to the fact the a discrete Majorana zero mode 
must be present inside the ring when the enclosed flux 
$ = ncf>Q (ip = 27m), when n is odd, but is absent when 
n is even. In order for the zero mode to appear or dis- 
appear, there must be a point where the gap vanishes 
and the transmission is perfect for 27m < ip < 2ix(n + 1). 
The perfect resonant transmission is thus a specific sig- 
nature for a gapless ID Majorana mode on the junction. 
However, the transmitted Majorana mode does not carry 
charge. We will see in the following that the transmitted 
Majorana mode leads to a step in the average current and 
a peak in the current noise. 



B. Current and Noise 

We now develop general formulas for the electrical cur- 
rent and noise in our geometry. Similar calculations have 
been performed previously in Refs. 17-29, We must com- 
bine the transmitted Majorana mode with an additional 
Majorana mode in each lead as well as properly treat the 
remaining incident electron channels. For each channel, 
the Dirac fermion electron operators may be expressed 
in terms of a pair of Majorana operators, 



a, a a, a . ■ a.a 

= 1e + tr lE > 



"-E 



~a,a . ■ ~a,a 

1e + lT l E > 



(2.7) 
(2.8) 



where c E a {c a E a ) describe modes incoming (outgoing) 
from lead a and channel a. We should note that, in gen- 
eral a lead a with A^ Q electron channels will have 2JV" 
electron and hole channels constrained by particle-hole 
symmetry. We are free to define our modes such that at 
= 7T, 7 lQ is the extension of the perfectly transmitted 
mode into the leads. Thus, there is 1 transmitting Ma- 
jorana channel and 2N a — 1 reflected Majorana channels 
with, in general, no additional constraints. 

We can express the relationship between incoming and 
outgoing Majorana modes in terms of a scattering ma- 
trix, 



which Se -> (t v Se& v , so that t. 



tE and 



Je = b E 1e > 



(2.9) 
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where a and b are now indexes for lead, channel, and 
Majorana type (7 or 77) and noting that S% b = due 
to particle-hole symmetry. Se has the general structure 



Se = 



i E 



(2.10) 



which allowing for only a single transmitting Majorana 
channel becomes: 



K. = 



i E e 



e > 



(2.11) 
(2.12) 

(2.13) 



where 
cients. 



and t E are single-channel scattering coeffi 
is a (2N a - 1) x (2N a - 1) dimensional Majo- 
rana reflection matrix representing the remaining chan- 
nels, and eft- = S ai Sbj. Our plots use for r' 



model scattering coefficients in Eq. (2.3) and (2.4 1. S 



E an d t% the 



and all other matrices for our model are 2 (A 1 + A 2 ) x 
2(iV 1 + A 2 ) dimensional. 

The operator for the current flowing out of contact a 
is given by 



~a,at ~a,a 
C -E C E 
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where A% 



S E T,ySE, Je nere i s an element of a 
column vector of Majorana operators in channel a, and 
= P a (£) a y , where P a is a projector onto the modes 
in lead a and a v is the Pauli matrix coupling 7°'" and 
r] a ' a . Omitted details for this calculation are presented in 
Appendix [X] The average current I a = (I a ) thus reads 
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(2.15) 



in which we have used the following definition 



Qe = (1-e1e) 
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(2.16) 



[f^P^ + f-E^y 

P 

where f E ' = f p {E) ± [1 - f p (-E)], the sum and differ- 
ence of the electron and hole Fermi functions. 

Similarly, the low-frequency noise power P a @ — 
j+™ dt((I a (t)i' 3 {0))-I a I^) can be calculated as follows 
using Wick's theorem. We find 



P 



L 

2e 2 v F 



E.E' 
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lE ,)Af a A^ 
J2^[AeQ-eA e Q t e }. (2.17) 
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FIG. 3: Conductance as a function of ip in the single-channel 
limit. The green curves represent the ideal case for zero tem- 
perature and zero voltage in each lead with Ae = 0.05Ao. The 
purple, orange, blue, and red curves respectively represent the 
cases for Ae/A = 0.5, 0.2, 0.1, and 0.05 with T = 0.005A 
and eVi = eV 2 = eV = 0.05 A . 



The detailed derivation is explained in Appendix [B| 



III. EXPERIMENTAL SIGNATURES 

Here we provide a description of the current and noise 
observables which characterize the ID gapless Majorana 
channel in the J J. We begin with the average current and 
conductance at each lead, finding that they are indepen- 
dent of the applied voltage at the other lead and display 
sharp steps in the low temperature and small voltage 
limit. We then consider the noise power across the leads 
and at the same lead. The cross noise signal contains 
only one term which exhibits a peak at (p ~ it, though 
is suppressed by interchannel scattering. The diagonal 
noise displays a more complicated signal but exhibits a 
peak that persists even in the large A 1 ' 2 limit. 



A. Average Current 

We will begin with a calculation of the average electric 
current in the limit that there is only a single channel 
in the electrical contacts. In principle, this could arise if 
there was a quantum point contact separating the leads 
from the surface states. Even away from this limit, how- 
ever, the simplicity of the result will aid the understand- 
ing of the more general results, which we present in the 
following section. 



1. Single- Channel Limit 

In the single channel limit, there is only one addi- 
tional Majorana mode r; 1 '" i n eac h lead, so that W — 
±1. Using Eq. \2.1ty and \2A7\ and setting Y,e = 
L/(2ttUve) J_°° dE, we obtain 



r = 



2h 



dE(l-r E )f«(E). 



(3.1) 
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In Fig. [3] we plot the conductance G — I /V as a func- 
tion of phase <p predicted by Eq. (3.4) for several values 
of Ae. G(ip) exhibits sharp steps at ip = 2ir(n + 1/2), 
provided Ae <C Aq. This requirement is equivalent to 
having the coherence length £ ~ Tiv E /Ao <C L, so that 
quasiparticle tunneling across the superconductor is sup- 
pressed for if ~ 0, leading to perfect normal or Andreev 
reflection, r(ip ~ 0) = ±1. In this limit the width of the 
step is determined by the maximum of Ae and T . For 
T <C Ae <C Aq the linear conductance is simply, 



G{ip) = (l-r Q )e 2 /h 



(3.2) 



where from Eq. (2.6) ro = tanh[(Ao/Ae) cos ip/2]. r$ 
switches between —1 and +1 over a range Sip ~ Ae/A <C 
1, so that over that range G exhibits a step 



AG = 2e 2 /h. 



(3.3) 



In the following section we will show that when there 
are additional channels, the step is still present, but its 
magnitude is suppressed. 

Eq. (3.2) can be understood in terms of the four 



elementary scattering processes for a particle at the 
Fermi energy incident from one of the leads. The 
probability for reflection as an electron (or hole) is 
|1 + (— )ro| 2 /4, while the probabilities for transmission 
as an electron or as a hole are both (t ) 2 /4. For 
Ae « T « A , G(ip) = G(m v /T)e 2 /h, with G(X) = 
l-Xjl dzVT^/ cosh 2 pfz). 



2. General N- Channel Current 

In the general case of many electron channels in each 
lead we find 

= ± £J dE f^i^N* -2Re[r%*f^ E ] 

-Tr[Z y ir E '%r E ]}. (3.4) 



This expression contains terms that do not depend on 
the scattering of the mode that is perfectly transmitted 
at ip = ir. In general, these terms will depend on ip, 
but they will not exhibit the singular ip dependence as- 
sociated with the critical mode. Thus, we extract the 
singular terms, which we denote by 7", that depend on 
r a and t a . 



I a = - — 
1 2h 



dE fl 



(3.5) 



In the limit T, V a < Ae, we have S E ~ S . Since S E = 
S*_ E , So is a real, orthogonal matrix. The conductance 
jump at ip ~ 7r due to the step in rg then becomes 



AG a 



h r22 '° • 



(3.6) 



In general, f 22 depend on the details of the in- 
terface between the Josephson junction and the electrical 
contacts, and can vary in both sign and magnitude. We 
will not attempt to compute it in detail here. Rather, we 
will note that f 22,0 is an element of a (27V Q -1) x (2iV Q -l) 
orthogonal matrix. Under the assumption that all such 
matrices are equally likely, the typical value will be 



22,0 



0(1/. 



(3.7) 



Thus, though the magnitude is suppressed, the conduc- 
tance still exhibits sharp jumps in the limit T <C SE <C 
Ao. Importantly, the conductances G 1 and G 2 from the 
inside and outside leads should exhibit jumps at the same 
magnetic flux. If observed, these conductance jumps rep- 
resent a clear signal of a quantum phase transition in the 
system, corresponding to the insertion or removal of a 
delocalized Majorana mode in the superconducting ring. 
The magnitudes and signs of these jumps are each char- 
acterized by f22 q. We will see that the same parameters 
characterizes the signature in the noise correlations and 
that the cross noise correlation peaks occur at the same 
flux as the conductance jumps. 



B. Noise Power 

Next we calculate the current-current correlations 
which contribute to zero- frequency noise power P a P . We 
find that the cross correlation P 12 is characterized be a 
peak at ip ~ ir due to a single term which corresponds to 
Majorana transmission across a ID gapless channel and 
if detected at half integer multiples of (f>o gives an unam- 
biguous signature of a quantum phase transition and the 
insertion of a delocalized Majorana mode into the ring. 
For a single channel, the peak has a universal height, 
while, for many channels it is suppressed by the same 
factors that led to the suppression of the conductance 
steps. The diagonal noise P aa will be presented in its 
single- channel and many-channel limits with relevant be- 
havior highlighted. Unlike the cross noise, the diagonal 
noise contains singular terms which remain as N 1,2 — > 00. 



1. Single- Channel Limit 



channel limit we find 



Taking Eq. (2.17) in the previously-discussed single- 

(3.8) 



PLZ = 4hJ „ d^fE^Etl, 



/"' = [/; / din r) 2 f l lr- 



- a-r E ) 2 f E A f E A + \t E \ 2 fHfp 2 ]. (3.9) 

P 21 and P 22 follow from interchanging superscripts 1 O 
2 in Eq. (3.8) and (3.9), respectively. 
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FIG. 4: Cross noise power P 12 as a function of ip in the single- 
channel limit. The green curves represent the ideal case for 
zero temperature and zero voltage in each lead with Ae = 
0.05A . The purple, orange, blue, and red curves respectively 
represent the cases for Ae/Ao = 0.5, 0.2, 0.1, and 0.05 with 
T = 0.005A and eVi = eV 2 = eV = 0.05A . 
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Fig. |4] shows the non-local noise correlation P 12 (ip) 
evaluated for eV\ = eV 2 = 10T, for representative val- 
ues of Ae. For Ae <C Ao, P 12 (ip) exhibits a peak near 

= 2ir(n + 1/2). For Ae < T the peak height is sup- 
pressed by a factor ~ exp(— irT/ Ae). Observation of the 
peak in the noise correlations requires eV, T < Ae. 



In the limit eV,T < Ae Eq. (3.8) and (3.9) reduce to 



P 11 = 



2h 



T{ (to)" 



~2T~ 



F 



eV~ 
~2T~ 



P L 



= k T ^ 2 



(3.10) 



eV+ 



- F 



eV~ 
~2T 



(3.11) 



V 1 ± V 2 
= to 



and 



and rg 



where F(X) = Xcoth(X), V ± 
made the assumption that <q = 
ro. This assumption is generally valid as most exper- 
imental systems in this geometry will be adiabatically 
connected to our model system in this low-energy limit. 
This leads to a striking behavior in the zero-temperature 
limit. For T < Vi, V 2 < Ae, we find that 



FIG. 5: Cross noise and diagonal noise at tp — n in the single 
channel limit. The green curves represent the ideal case for 
zero temperature and zero voltage in each lead. The purple, 
orange, blue, and red curves respectively represent the cases 
for T/Ao = 0.1, 0.05, 0.02, and 0.005 with eVi = 0.05A and 
Ae = O.lAo. 



The current flowing across the junction, / = I 1 
noise P~ = (6I~6I~) which goes as 



't 2 \V 1 -V 2 \ 



-I 2 has 



(3.15) 



The features in G((p) and P ab (ip) near (p = tt, which are 
predicted to occur over a width Sip ~ max(Ae, eV, T)/A 
constitute a signature for the Majorana fermion modes 
associated with the Josephson junction. They are present 
because over the range Sip, a Majorana zero mode is 
transferred from one end of the junction to the other. 



2. N-Channel Generalization of the Noise Power 
When there are multiple channels the cross correlation 
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2h 



2/t 



(to) 2 |K„ax|, 

(to) 2 V min sgn(V m 



(3.12) 
(3.13) 



where V max (V m in) is the voltage V a=1,2 that has the 
largest (smallest) absolute value. Obviously the diago- 
nal (off diagonal) noise correlations are sensitive to the 
maximum (minimum) of the voltages of the leads relative 
to the superconductor. In Fig. [5] we plot the noise corre- 
lations at the peak m = as a function of V2 for fixed V\ 
and representative temperatures. The fluctuation in the 
total current I + = I 1 + 1 2 flowing into the superconduc- 
tor is given by P + = (SI + SI + ) which has the following 
simple form: 



-t 2 \V 1 + V 2 \ 



(3.14) 
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In the limit, eV,T Ae, the cross noise still maintains 



the form of Eqn. (3.11 ) but is suppressed: 

o2 



2T \ 2T 



f 1 r 2 

u 22,0" 22,0- 



(3.17) 

For T <C Vi, V 2 <C Ae, the Fermi functions become step 
functions and the cross noise becomes 



p 3 V 
2h 



Sgn(V,nax) t (ip) f 22 . r 22 n ■ (3-18) 



7 



which still maintains the same V m i n , V max dependence as 
3.13 The height of this peak as a function of ip is related 
to the conductance jumps in Eq. (3.6| by the scattering 
parameters ir^o which can be independently measured 
using the average current at each lead. Thus, for many 
channels we expect the cro ss noise to be suppressed by a 
factor of order {l/V^N 2 ). 

The diagonal noise signal P 11 has a more complicated 
dependence on ip as well as the elements of and is fully 
derived in Appendix[Cj However, we find that in the large 
TV 1,2 limit, in which the 0(l/y/N a ) terms are completely 
suppressed, there remains a peak in the noise that gives 
a robust signature for the transmitted Majorana mode. 
In particular, we find 



Ah 



/+oo 
dE f+^fp 1 (2N 1 -\t, 
-OO 



+ f+ifp 2 



t E \ 



(3.19) 



In the T, V <C Ae limit the singular piece of this becomes 



-(t ) 2 T{F 



F 



eV~ 
~2T 



The voltage dependence of P^ anv ~ is further clarified in 



the T < Vi, V 2 < Ae limit: 



P. 



many,"/ 



2h 



(to) 2 



\V m 



w 1 



(3.21) 



which vanishes if | V^ 1 1 < \V 2 \. This behavior is quite 
distinctive and gives a clear indicator of a single gapless 
Majorana channel. There will be a peak in the current in 
lead one due to an applied voltage in lead 2, but no peak 
if a voltage is only applied to lead 1. In all of these cases, 
the diagonal noise still contains singular pieces in the 
large N 1,2 limit and therefore provides, of all the quan- 
tities discussed in this paper, perhaps the most robust 
signature of gapless Majorana modes. 



suppressed by 1/N. The diagonal noise includes a peak 
that is not suppressed for large N. 

We now briefly discuss some relevant issues for exper- 
imentally implementing our proposal. The number of 
channels of the leads is an important parameter for de- 
termining the size of the singular contributions. For a 
ring geometry, as in Fig. [T] this can be roughly estimated 
as N a ~ kpR a , where R a is the radius of the inner or 
outer edge of the ring. To minimize this, it is clearly 
desirable to control the Fermi energy of the topological 
insulator surface states, such that the Fermi energy is 
close to the Dirac point. In this case, N a ~ EpR a /Uvf, 
where vf is the velocity of the surface states. For Bi2Se3, 
Hvf ~ O.leV nm, so for R ~ 1^, and Ep ~ lOmeV, 
N ~ 10. 

A second key parameter in our theory is the level spac- 
ing Ae = Hvf/L. To observe sharp features in the cur- 
rent and noise at ip — tt we require Ae <C Aq, so that 
L is larger than the coherence length £ = Tlvf/A. This 
ensures that for cp ^ tt the transmission of quasiparticles 
across the superconductor is exponentially suppressed. 
Since the noise peak is suppressed for Ae < T observation 
of non-local noise correlations requires that L not be too 
large. There is ample room to satisfy these constraints 
experimentally. For example, in Ref. [14] devices with 
Ti/Al electrodes (A ~ 150/ieV) were studied. While vp 
is not known exactly for these devices, an upper bound 
is the velocity characterizing the Dirac surface states of 
Bi2Se3, Tivf ~ O.leV nm. This leads to £ < 0.6/im. For 
longer junctions, in which L > Lp = hvp/T or the inelas- 
tic length L; n , the noise correlations will be suppressed, 
but the step in the conductance as well as the diago- 
nal noise peak remain robust, provided the flux through 
the ring (and hence the phase ip) can be controlled with 
the applied magnetic field. For this, it is desirable to 
minimize the self-inductance C of the ring, so that the 
Josephson energy, Ej cos 2ttQ /</>o , which tends to quan- 
tize the flux, is dominated by $ 2 /2£. In Ref. Q31 de- 
vices with L ~ 1/im had critical current I c ~ 1/iA. Using 
Ej = hl c /e and C ~ AirRlogR/L (for a ring of radius R 
and thickness L), we find that this condition is satisfied 
for R < lOO^m. 



IV. CONCLUSION 

In this paper we have computed the electrical current 
and noise for a Josephson junction structure on the sur- 
face of a topological insulator that allows a clear signa- 
ture of the gapless Majorana mode, predicted at phase 
difference <p = tt. We predict that the average current 
exhibits sharp steps as a function of phase difference for 
a long junction at low temperature and voltage. The di- 
agonal and off diagonal noise correlations exhibit peaks 
at ip = tt. The amplitudes of the singular steps and peaks 
are predicted to be universal in the case where the electri- 
cal contacts couple via a single channel. For N channels, 
the singularities remain finite, but the current steps are 
reduced by l/y/N, while the cross noise correlation is 
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Appendix A: Average Current Calculation 

For the purpose of clarity, we will work out in detail 
our derivation of the Af-channel current and noise, be- 



ginning with our choice of a Majorana basis and working 
in this appendix up to the average current and its single- 
channel and many-channel limits. Appendixes [B] and [C] 
will expand upon this work up to the general expression 
for zero-frequency noise power and its forms in the inter- 
channel and same-channel cases. 

We begin with the choice of constructing Majorana 
operators 



Te 
Ve 



\«3 



--E 



(Al) 
(A2) 



where a represents a given channel and these linear com- 
binations have been chosen such that these new operators 
are +1 eigenstates of the particle-hole operator S = t x K 
in the electron-hole basis of 



n't 



ce 
,t 

"-B 



(A3) 



where we should note that j E is a Majorana of type 7 (77) 
for odd (even) a. There will always be an even number 
of modes since we have artificially doubled the electron 
channels in each lead as to handle both normal and An- 
drecv reflections off the superconductor. The alternating 
pattern of entries in 7 b allows us to summarize Eq. (A5| 
as follows 



Qe = (1-e1e) 



E Qz ab 

0=1,2 



1 r 

4 



(A8) 
(A9) 



The outgoing operators 7^ = S E b j E where Se is a 
2(N 1 + N 2 ) x 2(7V 1 + N 2 ) dimensional scattering matrix 
that obeys the property Se = S*_ E i which can be derived 
by Hermitian conjugating je and using j E = 7_b, and 
is a consequence of particle-hole symmetry. This allows 
us to write the current operator in a much more compact 
form 



These operators obey the additional property that 7!^ = 
7b, such at E — they obey the Majorana relation 
7q = 7o- We note that since our new operators are just 
linear combinations of electron creation and annihilation 
operators, they still have canonical anticommutation re- 
lations {^eiie'} — <^ ab ^E,E' and still obey Wick theo- 
rem when calculating higher order correlation functions. 
The two kinds of Majorana have valid contractions with 
themselves and between species, leading to four correla- 
tors that will be of use: 

<7^b7b<) = W-eVe') = \fp ab SE,E> , (A4) 



(Y-evU = -(V-ElE') = -Ji?> ab 8E,E> , (A5) 



where f E = f(E) ± [1 — /(— E)] and comes from substi- 
tuting the definitions of our Majorana operators into the 
above correlators and noting that {c a E c b E ,) = /eS^Se^' 
and {c a _ E c h l E ,) = (1 - f- E )S ab S E ,E', the Fermi distribu- 
tions for electrons and holes respectively. From here, we 
can introduce the current operator 



r 



ev F sr^ , at 

-£-2>-* 

E 



(A6) 



E" is the matrix for charge- weighted momentum through 
lead a in the electron-hole basis and we have expanded 
ip E to be 2(7V 1 + N 2 ) dimensional for our two-lead ge- 
ometry where N a is the number of channels in lead a. 
The —E for ift comes from the delta function in en- 
ergy that we get by time-averaging and summing over 
individual operator energies. E" = -P a E 2 where P a is 
the projection matrix into the subspace of lead a and 
E 2 = t(N) <E> (T z . Rotating this into our Majorana basis 



r = 



evF 



E^ s 



a,ab„,b 

v 7b 



J-eW 7B : 



(A7) 



I a = 



ev F 



E a Aa,ab b 
1-E A E 7B • 



A% = E£ - S E ^S E . 



(A10) 
(All) 



With all of these definitions in place, we can finally 
begin to take the expectation value of I : 

jl _ /fl\ _ el ' F S^/~,a „,b \ A l ^ b 



L Y.^-e1e)A 

E 

E Q C E A E ab 
E 

E Tr ^i 



evF 

~J7 



evp 

1 



El j 



(A12) 



which gives us a general average current of 



4L 



^{/-^Tr^-E^E^B] 

E 

+/+' 1 Tr[P 1 4E^B] 
+/ £ ' 2 Tr[E 2 y 4E^B] 
+f E -' 2 TrlP 2 S E Z 1 y S E }}. (A13) 

We can introduce into this a general S-matrix of the form 



r E t| 
i E r| 



(A14) 



in which unitarity restricts rgV| + tjj Itjj = 1. The 
average current now becomes 

^/-' 1 {2A 1 -Tr[E,4 t E,r £ ]} 



ri = 

4L 



fp 1 Tr[r^E,r B ] + f E 2 Tr^E^I] 



(A15) 
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where we have used E y to represent any square matrix 
of the form 1 ® cr y and E^ to represent any projection of 
a matrix of that form. This means that our Y, y matrices 
will possibly be different sizes; nevertheless they are care- 
fully ordered in such a manner that all matrix products 
are still valid. 

In our Majorana basis, only a single channel may be 
transmitting, which restricts the form of our S matrix 



n_a ja „1 1 

t E = t E e 
r R = r E e 



(A16) 
(A17) 



where r E and t E are scattering coefficients for the sin- 
gle transmitting channel and the elements of the matrix 
<^ij = SaiSbj. The scattering matrices for the remaining 
channels r E ^ r E are in general unequal as the leads 
differ in channel number and structure and are otherwise 
unrelated by additional symmetries. These matrices have 
the structure 



T 
R E 



(A18) 



where is the zero column vector and R E is an unde- 
termined (2N a — 1) x (2N a — 1) dimensional reflection 



matrix. Noting that t E ^T, y t E 

while taking ^2 E = L/(2irhvF) f_°° dE, we arrive at a 
final answer for the total average current at the first lead 

e f + °° 

I 1 = Ih J_ dE f E ' l {2N x -2Re[r 1 E *r 1 22iE ] 

-Tr[E y r^^ y r E }} (A19) 

of which the part containing singular Majorana behavior 
as ip goes through it is 



and E„e n E 



■12 



3 



e 

2h 



/+oo 
dE /^{Re^*^]}. (A20) 
-oo 



If we take the single- channel limit, i.e. the case of a 
quantum point contact, N 1 = N 2 = 1 and r E is reduced 
to P22 e = — ^22 e = 1 an d the average current becomes 



e 

2ft 



+ 00 



(A21) 



In the more realistic many channel limit, the elements 
of r E ~ 0(l/\/N a ) but with random phases. This means 
that they will, in general, not add coherently such that 
in the large N a limit, terms that contain r E will die off. 
In this limit the average current reads 



I 1 

many 



e 

2h 



/+oo 
dE f^N 1 
-oo 

and the <p dependence is suppressed. 



(A22) 



Appendix B: Cross Noise Calculation 

In this section, we calculate in terms of our previous 
matrices the general expression for noise power, special- 
izing at the end to the case of cross noise. Zero-frequency 
noise power, P a P = /+~ dt({I a (t)P{Q))-I a P), can be 
calculated using Wick theorem as follows: 

pa0 = ^^{l-ElEl-E^ATA^ 
E,E' 

2 

e V F \ * s^ad/^bc Aab.a Acdj3 f^acr^bd Aab,a Acd,0 

— t / . We W-e a e a e ~ ^e h~e a e a -e 



L 

E 
E 



(Bl) 



From here, we can specialize to the a ^ j3 case and 
calculate the cross noise power. For a = 1, f3 = 2, 

2 

P 12 = 4^ ^[Se^SeQIeSe^SeQY] ■ (B2) 

E 

After specializing the S-matrix elements and taking 
Y, E = (L/(2nhv F )) J dE, P 12 becomes 



pl2 

7 



e 

ih 



/+oo 
dE f E ^f E ' 2 t E t 2 E r^ E r 2 ^ E , (B3) 
-oo 



where we have used P^ 2 to note that all of P 12 behaves 
singularly as the system goes through its critical point. 
In the single channel limit discussed in Appendix [A] 

PU = fft J „ dE f * (B4) 
In the many-channel limit, P 12 -> as O^/y/^N 2 ). 



Appendix C: Diagonal Noise Calculation 

The diagonal noise calculation follows very closely the 
cross noise calculation. However, unlike P 12 , P 11 does 



not have a simple general form. We begin with Eq. (Bl ) 



L ^ 

E 



Tv[A E Q_ E A E Q T E ], (CI) 



which given Eq. ( A17 1 and (A18) becomes 
„l r+°° 



e 
8ft 



dE 



/-i/s' 1 ! 4 ^ 1 - 2TV[£^IVi] - 2\t E \ 2 } 

f E ' l f E ' l { - 27V 1 + 2Tr[E„r^E„^] 
Tr Pi, r 1 ^ E y r E E a r^ f E y r E ] } 

Ue-'Ie' 2 - / B ' 1 /^ 4 ){2Tr[E,/^^E y ^E !/ ^]} 
f^fp 2 2\t E \ 2 ) (C2) 
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under the usual substitution for ^2 E . Further simplifying 
this, we can obtain expressions for the traces of scattering 
matrices: 

Tr^t^E^E^] = \t E \ 2 {r E E y ^) 22 
Trp^E^] = 2Re[r|?f 2 - 2iE ] + Trp^E^] 
Tr[E } ,ir^ t E y r^E l ,r^ t E 1 ,r^] = 

Tr[E^E a f^r|;E y f^].(C3) 



2 Re 



This very complicated expression does not have, like 
the current and cross noise, clearly separable singular 
pieces in its general N-channel form. However in the 
single-channel limit as described in Appendix [S] it sim- 
plifies significantly: 



Pn = 



4/; 



/+oo 
dE[\l-r\\ 2 f^f, 
-co 



.1 

E 



where we have exploited in the second term that f_g 



ft 



and that 



In the many channel limit, 



the diagonal noise becomes quite simple: 



P. 



li 



4/; 



+00 



dEf + if^{2N'-\t E \ 2 } 



+f^ft' 2 \t E \ 2 ,(C5) 



(1 - ?i) 2 / B 'V £ ' 1 + \tE\ 2 fHfp 2 } , (C4) channel limit if V 1 = V 2 



where unlike with P 12 , the <p dependence is mostly pre- 
served and is much more clearly extracted than in the 
general, N-channel case. Additionally, the singular part 
of the diagonal noise, P^ 1 , goes to zero in the many- 
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